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Abstract 

The phenomenon of atomic population trapping in the Jaynes-Cummings 
Model is analysed from a dressed-state point of view. A general condition 
for the occurrence of partial or total trapping from an arbitrary, pure initial 
atom-field state is obtained in the form of a bound to the variation of the 
atomic inversion. More generally, it is found that in the presence of initial 
atomic or atom-field coherence the population dynamics is governed not by the 
field's initial photon distribution, but by a 'weighted dressedness' distribution 
characterising the joint atom-field state. In particular, individual revivals in 
the inversion can be analytically described to good approximation in terms of 
that distribution, even in the limit of large population trapping. This result is 
obtained through a generalisation of the Poisson Summation Formula method 
for analytical description of revivals developed by Fleischhauer and Schleich 
[Phys. Rev. A 47, 4258 (1993)]. 
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I. INTRODUCTION 



One of the most fundamental models of quantum optical resonance is the interaction of 
a single two-level atom with a single quantised mode of radiation, described by the Jaynes- 
Cummings Hamiltonian Despite being simple enough to be analytically soluble in the 
rotating-wave approximation, this model has been a long-lasting source of insight into the 
nuances of the interaction between light and matter. It has led to nontrivial predictions, 
such as the existence of 'collapses' and 'revivals' in the atomic excitation 0], and has also 
allowed a deeper understanding of the dynamical entangling and disentangling of the atom- 
field system in the course of time Further interest in the Jaynes-Cummings model 
(JCM) comes from the fact that these predictions are directly accessible to experimental 
verification. A JCM interaction can be experimentally realized in cavity-QED setups |^ 
and also, as an effective interaction, in laser-cooled trapped ions (in which case the 
ionic harmonic motion assumes the role normally played by the field mode). For example, 
revivals in the atomic excitation have recently been observed in a cavity-QED experiment 
0, providing direct evidence for the discreteness of photons. 

In spite of these successes, a closed analytical description of the collapse-revival pattern 
has so far proved to be elusive; however, an elegant approximation scheme valid for a number 
of initial conditions has been presented by Fleischhauer and Schleich |Q, improving the 
earlier work of Eberly and co-workers 0. Among other things, they have demonstrated 
that, when the atom is initially completely excited or de-excited, and the initial photon- 
number distribution of the field (P„) is sufficiently smooth (as is the case in a coherent 
state) then the shape of each revival is a direct reflection of the shape of Pn- This direct 
relationship can be affected by the presence of initial atomic coherence. It has been noticed 
0] that if the atom is initially prepared in a coherent superposition of its energy eigenstates, 
then the revivals can be largely suppressed, effectively freezing the value of the atomic state 
populations. Appropriately, these initial atomic states have been referred to as "atomic 
trapping states" P,p!0| . This population trapping has been connected to the existence of 
a specific phase difference between the relative phase of the atomic superposition and the 
phase of the initial field state . In particular, in the resonant case exact trapping has been 
shown to exist [|ll| if the field is initially prepared in a phase coherent state [0] (eigenstate of 
the Susskind-Glogower phase operator [T^). As has been pointed out by Cirac and Sanchez- 



Soto this can be understood by the fact that in these cases only one state out of each 
pair of dressed eigenstates of the JCM is ever populated. 

We now note that, since the atomic excitation dynamics can be modified merely by 
altering the initial atomic state, then it is apparent that the initial photon distribution cannot 
be the sole responsible for the existence and shape of revivals. In fact, as we demonstrate 
in this paper, the key to understand the collapse-revival pattern under more general initial 
conditions is to consider not only the initial coherence of the atom and field by themselves, 
but also their joint properties as a single quantum system. This is so even if the system is 
initially disentangled. Building on Cirac and Sanchez-Soto's observation, we find that the 
revival structure depends essentially on the relative weight of each dressed eigenstate in the 
initial atom-field state. We are able to identify an atom-field variable which plays, in the 
general case, the same role as that of the photon distribution when the atom is initially 
excited or de-excited. This allows us to estimate for any initial condition the degree of 
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suppression of the revivals (that is, the amount of trapping which occurs). In addition, 
in some particular cases we are able to describe analytically to a good approximation the 
shape of the revivals, even when they are partially suppressed. To illustrate the dressed- 
state distribution dependence of the revival patterns, we show that essentially identical 
examples can be found both in the presence of initial atomic coherence, and in cases where 
this coherence is not only atomic, but where atom and field are already entangled. 

Our work is organised as follows: in section II, we rewrite the dynamics of the Jaynes- 
Cummings model from a point of view based on the entangled dressed-state eigenbasis. With 
the help of an appropriate coordinate system, the expression for the atomic inversion assumes 
a simpler form; as a consequence, we are able to calculate bounds for the evolution of this 
quantity for arbitrary initial conditions. In order to illustrate our technique, specific cases 
are treated in section III. In section IV, we generalise Fleischhauer and Schleich's scheme to 
include the cases where trapping occurs. Our conclusions are summarised in section V. 



II. DRESSED-STATE COORDINATES FOR THE JCM 

A. The Model 

The Jaynes-Cummings Hamiltonian, on resonance and in the rotating-wave approxima- 
tion is given by 

H = -hcoaz + hua)a + hX{a+a + a^(T_) (1) 

where u is the atomic transition frequency, A is the atom-field dipole coupling constant, a^(a) 
is the field photon creation (annihilation) operator and a^, a± are the atomic inversion, rising 
and lowering operators. The first two terms in this expression describe the internal energy 
levels of the spin-l/2-like 2- level atom and the harmonic oscillator- like field mode, and the 
interaction term can be straightforwardly understood as a one-excitation exchange between 
these two systems. 

An important feature of this fully-quantised matter-radiation interaction is that its 
steady states (known as 'dressed states') are entangled. Switching to an interaction-picture 
representation for convenience, it can be shown that these states and their corresponding 
energy levels are [T^: 

\n±) =l^[\e,n)±\g,n+l)] (2a) 



V2 



En± = ±hXVn + 1 = ±-n„ (2b) 

{Qn is the Rabi frequency). It is apparent that in all these cases the atomic state is completely 
undetermined, being an equal mixture of the excited and ground states. Thus, if we assume 
that at f = the atom and field are independently prepared in the states 

IV'(O))^ =p\e)+q\g) (3) 

oo 

l^mp =Ecn|n) (4) 

ra=0 
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then the initial state (0))^^ = |-0 (0))^ ® |$ (0))^, will evolve into 
\'^{t))AF =qco\g,0) + 

pCn COS (^Xty/r 



+ E 

n=0 

oo 

+ E 

n=0 



'n + 1 
—iqcn+i sin (At-v/rT+T 



gc„+i cos (^At^Ai+T 
—ipcn sin (^AtV^ + l) 



|e, n) + 



\g,n + l) 



(5) 



Despite being straightforwardly solvable in this way, the JCM is well-known for the fact 
that the time-evolution of most expectation values is usually expressible only in series form. 
For instance, assume that at time t = we have: 



1^(0)) = {p\e) + q\g))^(Y.Cn\n)] 

\n=0 I 

oo 

= gco 1^, 0) + ^ pCn |e, n) + gc„+i l^f, n + 1) 



(6) 



n=0 



The time evolution of the atomic population inversion is then given by: 



oo 

2 ^ \pCn COS {\t\/n + — iqCn+\ sin (^AtV^ + l) 

n=0 

oo 

- kcol^ + ^ (bc„|^ - |gc„+i|^) cos (2AtV?2 + l) + 



n=0 



-F2/m ((pc„)* gc„+i) sin (2AtV' 



'n 



where 



El |2 I I |2 1 



I^Col 



(7) 



(8) 



n=0 



In particular, if at f = the atom is completely excited {q = 0) or in the ground state 
{p = 0), and admitting for simplicity that Cq is negligible, then this expression reduces to 
the simpler form 



oo 

(a,) (t) = ±Y, Pn COS (2AvArTTt 



n=0 



(9) 



where F„ = |c„|^is the initial photon distribution of the field. Due to the y/n scaling of 
the various Rabi frequencies in this sum, no closed form for it is known. Nevertheless, it 
is apparent that features of this time evolution, such as collapse-revival phenomena, must 
be directly related to the characteristics of P„. The Fleischhauer-Schleich approximation 
scheme, which we shall discuss in section IV, allows this relationship to be explicitly described 
in certain situations, such as when P„ is sufficiently smooth. 

Apart from the simple case represented above, in general there is no way to relate the 
inversion dynamics exclusively to the initial field state. As can be seen from expression (^, 
it will usually depend on atomic and field variables in an apparently complicated way. In 
the next section, we show how this evolution can be recast in a simple form, reminiscent of 
expression (H), independently of the initial conditions. 
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B. Dressed-State Coordinates 



Our discussion will be based on the relative contribution of each dressed state for the 
global atom-field state. For simplicity, we restrict ourselves to pure initial conditions. There- 
fore, we may expand an arbitrary initial state in terms of |^±) : 



where 



n=0 



COS j |n+) +e-*^" 



sin ( y 1 1?^-) 



(10) 



(11) 



with the normalisation condition 10"^ = !. The parameters tf„ G [0, 1] , ^„ G [0, tt] and 

Xn, 4>n ^ [0, 27r] , will henceforth be referred to as dressed- state coordinates. 

The use of this coordinate system allows us to obtain a simple geometrical picture of the 
manifold of pure atom-field states. For each value of n > one can imagine a Bloch-type 
spherical shell corresponding to the two-level system formed by \n^) (Fig. 1), which are 
parametrised by the spherical coordinates 9n and (pn- The poles (6'„ = 0,7r) are associated 
to the dressed states |n±), and points on the equator (9^, 



to the states of the form 

\n+) + e~^^" |n— ) , which include the particular product states |e, n) and \g, n + 1) (cor- 
responding to 0„ = 0, 7r). Therefore the coordinate 6n gives a measure of what we shall call 
the ^ dressedness^ (or degree of proximity to the nearest dressed state) of the components 
\ipn) of |\E'(0)). We note that, while this property is related to the entanglement of {ipn) 
(since the dressed states are in fact maximally entangled), they are two different concepts: 
states with the same dressedness, such as |e, n) and cos /i |e, n) — i sin fi\g,n+ 1) can have 
different amounts of entanglement (as measured by the von Neumann entropy of their sub- 
systems). Meanwhile, the weight factors Wn measure the relative importance of each of these 
components. In physical terms, their squares w"^ correspond to the probability distribution 
for measurements of the total excitation number operator 



Xtot = + N. 



(12) 



on the initial state (0)) { corresponds to the probability for n -|- 1 excitations). We 
note that, in the special case where the atom is prepared in one of the states |e) or l^f), then 



9ri 



and 



(tt) for all n, while reduces to the photon number distribution: 



P^ = \{n\^m^f (if |e)) 



71+1 



|(n + l|<|.(0))^r (if \9)) 



(13) 



We now show that the dressed-state coordinate system allows a simple description of 
the time evolution of any initial state. From equations (^,b), the evolution of state ([To|) is 
given by 



\9,0) + 



00 

n=0 



0„ 



COS 




(14) 
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We can re-express this directly in terms of the dressed-state coordinates: 

On (t) = en (0) (15) 
Xn it) = Xn (0) - -^nt 
4>n if) = 0„ (0) - Vtnt 

Thus, Wn and On are both constants of motion, while the angles Xn and 0„ precess with 
a constant angular velocity proportional to the Rabi frequency A way of visualising 
this is by means of a geometrical image as in Fig. 1. The simplicity of this picture reflects 
the symmetries of the resonant JCM: the conservation of Xtot and the fact that, in the 
interaction picture, the energies in each 2-dimensional eigensubspace generated by |t?-±) are 
always symmetric with respect to zero. 



1. Population Inversion in Dressed- State Coordinates 



Let us now apply this coordinate system to describe the evolution of atomic variables. 
In the basis {|e) , l^f)} , the reduced atomic density operator can be expressed as: 



PA {t) 



Pee (t) Peg (t) 

iPlgif) Pggif) 



where 

Pee (t) 
Peg (t) 



1 oo 

- (l - + sin (On) cos (t)) ; Pgg (t) = 1 - Pee (t) 



n=0 



V2 

oo 

2 



COS 



+ e 



-i4>o{i) 



-AX-a+l{t)~Xn{t)) 



+ - WnWn+ie'"^"^''^' I COS 



ra=0 



X COS 







sin ^ 








V 2 / 



+ 



'ft, 



n+l 



- e^-^-W sin 



(16) 



(17a) 



(17b) 



We thus see that the evolution of the atomic level populations, and therefore of the 
population inversion, can be expressed in a simple form for any initial condition: 



(a, {t)) = 2pee (t) - 1 



-w 



-1 + H sin (On) COS (0„ (t)) . 
n=0 



(18) 



This should be contrasted with the much more complicated expression (|^) for this same 
quantity, written in terms of separate atomic and field coordinates. In fact, we can see that 
the present expression is a direct generalisation of eq. (^, obtained by phase-shifting the 
^th iqj-j^ jj^ f^j^Q sum by 0„ (0) and by replacing the photon number distribution P„ with 



sin 



(19) 
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It is not difficult to understand the reason for this substitution. First of all, as was noted 



by Cirac and Sanchez-Soto [|14|, dressed states have no population inversion. Therefore, 
if a component \tl)n) of a given state (0)) has a large 'dressedness' (sin (6'„) — >• 0), it 
will contribute very little to the overall inversion. Moreover, even if \il}n) does have a large 
inversion, this may still be of little consequence to the total average {az if)) if its relative 
importance w"^ is small. Thus, the product D„, which we shall refer to as the 'weighted 
dressedness distribution' gives the appropriate magnitude of the contribution of \il)n) to the 
inversion. In other words, for general initial conditions it is this distribution, not P„, that 
will control the evolution of the inversion. For instance, we can expect the existence of 
collapse-revival structures if D„ is 'narrow', in the sense that few values of n (and therefore 
few Rabi frequencies feature significantly in the sum. 



2. Population Trapping 

The main difference between the general evolutions allowed by equation ([TB|) and those 
obtained in the special case (^ lies in the fact that the distribution P„ is normalised, while 
Dn is generally not. This implies that, in general, there will be an upper bound for the 
allowed range of variation of the population inversion, given by: 

oo 

az{t))+w\\<Y.Dn = M <l (20) 

n=0 

As we can see, this bound is specified essentially by the average dressedness of the com- 
ponents of (0)). For example: if the most important components (those with largest 
weights w'^ are highly dressed (sin (6'„) is small), then the amplitude of variation of the 
population inversion will also be small (M ~ 0). In other words, these are the conditions 
for the existence of population trapping in the (resonant) JCM. The steady-state value of 
the inversion is given simply by — w^i, i.e., by the fixed population in the ground state. We 
emphasise that, given any initial condition, the bound M can be immediately calculated 
from the constants of the motion m„, thus giving an estimate of the amount of trapping 
that can be expected from that state's evolution IT 



III. EXAMPLES 

In order to illustrate our technique, we analyse three classes of initial atom-field states 
for which population trapping occurs: 

A. Perfect Trapping States 

In iTTH , Cirac and Sanchez-Soto found a class of factorised initial conditions which exhibit 



perfect trapping, that is, for which the atomic populations are strictly constant over time. 
These were of the form: 

I 1 °° 

I* (0)) = ^l-\z^\{z |e) + \g)) ® ^" 1^) (21) 

1 + n=0 
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where \z\ < 1. In this case, the field is prepared in an eigenstate of the Susskind-Glogower 
phase operator V = Z]^ol^)(^+l| EH ('phase-coherent state' fl^), while the atomic 



superposition is chosen so that the phase of the atomic dipole matches that of the eigenvalue 



z of the field state ||T9[. The resulting population trapping can thus be attributed to a suitable 
matching of atomic and field parameters. However, as was later noticed by the same authors 



1^ , a clearer explanation can be found by using a dressed-state point of view: Rewriting 



this state in terms of the dressed-state basis, we find 



1^(0)) 



1 - 


z^ 




1 + 


z^ 





\g,0) + Y.V2z''\n+) 



n=0 



(22) 



The reason for the existence of trapping now becomes immediately clear: all the components 
\ipn) of this state are completely dressed (sin (^„) = 0), and thus do not contribute to the 
inversion. 

Cirac and Sanchez-Soto also stated that these are the only factorised states exhibiting 
perfect trapping. We note that this is not strictly true: the same result is still achieved for 
all states of the form 



1 °o 

I* (0)) = ^l-\z'\iz \e) + \g)) ® E J H 1^) 

'1 + \Z^\ n=0 



(23) 



where j (n) can be ±1 and < 1 (in this case, the field state is not in general an eigenstate 
of V"). Once again, this can be seen by rewriting the state in the dressed-state basis. In this 
case, we obtain 



1^(0)) 



1 - 


Z^ 




1 + 


z^ 





n=0 



(24) 



where for each n there is a single dressed state \n±) present, whose sign is the same as that 
of the ratio ^^f^p-- The converse statement is also easy to show: suppose I^P (0)) is a perfect 
trapping state, of the form: 

oo 

1^(0)) =k{z\e) + \g))0Y.^n\n) 



n=0 



kao \g, 0) + E \e,n) + a„+i \g, n + 1) 

n=0 



(25) 



where /c is a constant and arbitrary. Then for each component of this state 

to be perfectly dressed it is necessary that: 



and hence 



a„+i = ±zan 



ttn OC ±Z'^. 



(26) 



(27) 



We recover expression (^) by noting that \z\ must be less than 1 for (0)) to be nor- 
malisable. Thus, for the resonant JCM there can be no population trapping with positive 
population inversion. 
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B. Zaheer-Zubairy 'atomic trapping states' 



While the initial conditions considered above lead to perfect trapping, preparing the 
appropriate field states requires elaborate state-engineering techniques. It has been shown 
(numerically) by Zaheer and Zubairy, however, that approximate population trapping is also 
possible if the cavity field is initially in a coherent state which are easily accessible in 
experiment This is portrayed in Fig. 2: if the atom is initially completely excited, the 
inversion evolves according to the familiar collapse-revival pattern 0]; however, by suitably 
rotating the atomic state, the revival peaks become more and more reduced, practically 
disappearing when the atom is prepared in state 



i^mA=^{\e) + e-^\9) 



(28) 



(where is the phase of the field's coherent amplitude a). Close inspection also reveals 
that, in this limit, the revival peaks become split into doublets. 

In their work, Zaheer and Zubairy gave no quantitative explanation for this quenching 
of the revivals, attributing it qualitatively to a "destructive interference between the atomic 
dipole and the cavity eigenmode ". Later, Cirac and Sanchez-Soto pointed out that the 
existence of approximate trapping in this case should be expected due to similarities between 



the properties of coherent states and of the phase coherent field states [|TT . 

With the help of our dressed-state coordinate system, we can now give a quantitative 
description of this effect. Let us suppose that the atom-field system is prepared in the state: 



^zz(tt,7,0)= COS7 |e) + e *^sin7 l^f) ® \a) 



(29) 



where the atom is in a coherent superposition of ground and excited states and the field is in 
a coherent state (a = |a|e*'^"). Expanding this state in terms of the dressed-state basis, we 
find that its dressed-state coordinates have the following values (we have only listed those 
relevant for the atomic inversion): 



sin On (a, 7)0 
cos (f)n (a, 7,0 

sin 0„ (a, 7,0 



[n + 1) cos^(7) + |ap sin^(7) 



{n + l)\ 

(n + 1) cos2(7) - |a|2e*2K-«) sm^{-f) 



(n -|- 1) cos^(7) -|- 
{n + 1) cos^(7) — 


a 
a 


^ sin^(7) 
^sin^(7) 


{n + 


1) cos^(7) — 
^{n + 1) sin 


^^2Qi2{ua 5)sin^(7) 
(27) sin {ua - 



(n -|- 1) cos^(7) — |Q;Pe*2('^'^~^) sin^(7) 



n> -1 



(On e [0,7r]) 



(30) 



e [0,27r]) 



The weight factors w"^ follow a Poisson distribution when 7 is equal to zero (they reduce to 
the photon number distribution of the coherent field) and have small deviations from that 
distribution for any value of 7. For moderately large therefore, w"^ is maximised around 

1 2 1 2 

the integer n^ax closest to |q;| or |q;| — 1. Meanwhile, it is not difficult to show that the 
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sin 9n distribution has an absolute minimum, located around rimin ~ |<^| tan^7 — 1, with 
minimum value given by 

sin(^„.„,J^sin(|z/,-e|) (31) 

This indicates the component of |\E' {a, 7, ^)) having the largest dressedness. 

Following the procedure outlined in section [11 B 2|, we can now see that to attain the 



highest possible degree of population trapping we must maximise the largest dressedness, 
and match it with the largest weight factor: 

rimax = rimin ^ tan7 ~ 1 (32a) 
sin(^„„,J = 0^1^^ = ^ (32b) 

We thus explain the "atomic trapping state" (p8|) found by Zaheer and Zubairy: it is the one 
for which the total atom-field state (a, 7, C)) most closely resembles a "perfect trapping 
state", in the sense that its most important components l^/jn) closely resemble dressed states. 

It is clear from expression (^) that the degree of trapping experienced from initial state 
l^zz (a, 7, 0) depends both on the relative weight 7 of each atomic state and on the phase 
difference (z^q — ^) between the atomic dipole phase and the field's coherent amplitude. In 
Fig. 3 we illustrate the change in the weighted dressedness distribution D„ as a function 
of these parameters: on the left-hand side, we take equal weights for |e) and \g) , satisfying 



condition (|32a|) , and vary (i^a —■ ^om | to 0. We can see the quenching of the distribution 
as the dressedness of the component |'0nmax) is increased. On the right-hand side, we keep 
the phase difference null and vary the relative weight from to |, obtaining a similar effect 
(in both cases, we have chosen a = 7). Comparing this figure with Fig. 2, it is apparent 
that there is a striking similarity between the profile of the D„ distribution and the envelope 
of the revivals in the corresponding time evolution, including the appearance of a doublet 
structure in the limit of the optimal trapping conditions. As we shall see in section IV, this 
is no coincidence. 



C. 'Even-odd' entangled states 

The achievement of population trapping in the previous examples was seen to be ulti- 
mately caused by the quenching of the D„ distribution, and therefore by the joint properties 
of the atom-field system. However, since the initial states considered were factorised into 
atomic and field states, it could also be argued that the trapping was due to a suitable 
matching of independent atomic and field parameters (i.e., the phase of the field's coherent 
amplitude matching the atomic dipole phase in the previous example). In order to stress 
the underlying importance of the dressed-state point of view, we now show that essentially 
identical evolution and trapping patterns for the population inversion can be obtained even 
when the atom and field cannot be considered as independent quantum systems. 

Consider the following class of states of the atom-field system: 

I'^EO (tt, 7, 0) = cos(7) |e)^ \even) ^ + sin(7)e*^ \g) \odd) ^ , (33) 

where \even) ^ = (|«) + \~ct)) and \odd)j: = -^(|tt) — |~Q^)) are respectively even and 
odd coherent states of the cavity field (here |a| is assumed large enough so that these states 
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can be considered normalised). It can be seen that, apart from the hmiting cases where 
7 — s> 0, |, these states are entangled. Furthermore, since \even) j: (^\odd)jr^ have nonzero 
amphtudes only for even (odd) photon numbers, an expansion of this state in terms of 
dressed states |n±) features only those with even n (in other words, W2n-i = 0). Physically, 
this implies that states of this form have zero average electric and magnetic fields, and zero 
average atomic polarisation. 

(^a + d^) = (a - a"^) = (a^.) = (ay) = (34) 

Thus, if at t = the atom-filed system is in state {"^eo {<^,l,0)y ^^e population inversion 
evolves according to eq. (plSf), but with only the even terms present: 



{^z (t)) = ^In sin(6'2„) cos(02n(^))- (35) 

n=0 

Now, it turns out that all these even-indexed dressed-state coordinates have values vir- 
tually identical to the ones listed in expression (^) above for the state \^zz (a, 7,0) with 
the same parameters a, 7, ^. The only difference is a multiplication of the weights W2n in the 
present case by a factor of 2, to preserve normalisation in spite of the absence of odd- indexed 
components. The resulting time-evolution of the atomic inversion, plotted in Fig. 4, shows 
that the quenching of revivals and appearance of a doublet structure as {fa — i) ^ ^ pro- 
ceed in essentially the same manner as in the previous case, the main difference being the 
doubling of the frequency of revivals [|20|. We note, however, that in the present case neither 



of the parameters Va,^ can be unambiguously assigned exclusively to the atom or field, so 
that population trapping in this case must necessarily be understood as a result of the joint 
atom-field properties of the state in question. 

It is also worth remarking that, due to the lack of atomic polarisation, the atomic reduced 
density operators of these states are always diagonal in the |e) , \g) basis . Therefore, the 
reduced entropy Sa = —Trpa In pa, which measures the degree of entanglement between the 
atom and the field, is entirely determined by the population inversion: 

Sa (t) = -Pee (t) In (pee (t)) - (1 - Pee (t)) lu (1 - p^e (t)) 

= -l{l-{aAtmn Q(l-(a.(t)))) 

-^(l + (^.W))lnQ(l + (^.(t)))). (36) 

Thus, in this case the existence of population trapping is equivalent to the atom and field 
remaining (nearly) maximally entangled during their entire time-evolution. For any given 
state l^^EO («5 1, 0) 5 ^ lower bound to the value of Sa (t) at any given instant of the evolution 
is then given by: 

^min = (1 - M) In (1 - M)) - 1 (1 + M) In (1 + M)) (37) 

where M is the bound defined in eq. (^). For instance, in the case of the 'trapping state' 
depicted in Fig. 4, where, a = 7, 7 = ^,(^ = 0,we have Smin = 0.69005, very close to the 
maximum possible value ln2 ~ 0.69315. Finally, we note that despite the existence of this 
lower bound, it is still possible to devise schemes by which such entangled atom-field states 
can be constructed [ITBllTTi. 
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IV. POISSON SUMMATION FORMULA FOR REVIVALS IN THE CASE OF 

POPULATION TRAPPING 



In 0, Fleischhauer and Schleich obtained approximate analytical expressions for the 
evolution of the atomic inversion, using a stationary-phase method based on the Poisson 
Summation Formula [^. They showed that in many cases the inversion can be written as 



a sum in which each term (t) is non-negligible only during a certain extension of time. 
This is in contrast for instance with the expressions used in section [11 A] above, where each 
term in the summation is periodic. 

Their work was restricted to initial states of the form \g)j^ 1'^')^', in which the inversion 
is given by eq. @. They were able to show that, if the photon distribution F„ of \^)p is 
sufficiently smooth, then the k^^ term in the Poisson sum for the inversion is given by 

„/ \ (XH- n\ 

= ^("=i^ -h^I-Tl (38) 



where P (n) is a continuous interpolation of P„. 

The main interest of this formulation is the fact that, if P„ is also sufficiently narrow, 
then the term Uk {t) describes to a high accuracy the evolution of the inversion during 
the k^^ revival. This can be readily seen from the formula above: this term describes a 



rapid oscillation in time, modulated by an envelope centred around t ~ ^""^v ^^^^ whose 
format is essentially that of P„. Therefore, for narrow enough P„, some of the terms Uk (t) 
become completely disjoint from the rest, describing an independent revival. Fleischhauer 
and Schleich were able to use this formulation to derive many interesting results, such as the 
decrease in amplitude and increase in width of each successive revival, and also the number 
of revivals that can be resolved for a given state, before they become scrambled due to their 
increasing width and consequent interference with each other. Finally, they were also able 
to extend the technique to some cases where P„ is not smooth, such as in squeezed states. 

Using our dressed-state formulation for the population inversion, eq. ([T8|), we are able to 
extend the Poisson Summation Formula method to even more general initial states, including 
ones with atomic or atom- field coherence. In the Appendix we show that, under suitable 
conditions, the behaviour of the atomic inversion during its k^^ revival is described by: 



Xt 



D{n) cosU„ (0) + ^ 



(39) 



where D (n) is a continuous interpolation of the 'weighted dressedness' distribution of 
the initial state. Thus, in general it is this distribution, not P„, that is reflected in the 
shape of each revival. In particular, in the case of population trapping, the quenching of the 
revivals mirrors that of D„. This explains the similarity between the proflles in Figs. 2 and 
3, including the doublet structure of the revivals in the limit of maximum trapping. 

In order to display the accurateness of this expression, we use it to calculate the revivals 
in the case of the Zaheer-Zubairy states \^ zz («, 7, 0) introduced above. From expressions 
(^) for the dressed coordinates, we see that in this case the weighted dressedness distribution 
is given by 
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Dn {a, 7, = V Qi in) + Ql {n) - IQx in) ■ Q2 (n) ■ cos 2 {v^ - i) 



(40) 



where 



Qi (n) = exp (- |a|^) cos^ (7) 

|2(n+l) 



l2n 



(^2 {n) = exp 



(n+l) 



sin (7) 



(41a) 
(41b) 



For sufficiently large a, these Poissonian distributions are well approximated by Gaussians: 

cos^ (7) 



Qi (n) 



2n \a\ 



exp 



sin2(7) 
Q2 (n) ^ , exp 



271 \a\ 




(42a) 
(42b) 



Extending these to continuous values of n and substituting in (^) (|39D we obtain an an- 
alytical expression for the fc*^ revival in the inversion. In Fig. 5 we plot the first two of 
these [k = 1,2) in the cases of maximum and minimum population trapping, alongside the 
corresponding exact evolutions obtained numerically. Despite a little distortion, agreement 
is seen to be very good. 

A similar calculation is also possible in the case of the entangled 'even-odd' states 
|^Eo(a,7,0) presented in the previous section. In this case, formula (|39| ) given above 
for the inversion is not valid, due to the strong oscillations in D„ [D^ = for odd n]. 
However, it is possible to adapt our calculations to this situation (see section |A 1| in the 
Appendix), obtaining the expression: 



it)) 



Xt 



,7r 



D (m) cos 



hmiO) + + Trfc - 



2nk 



4 



m+l: 



_ A2t2 



(43) 



for the inversion during the A;*'' revival. Here D (m) is a continuous interpolation of the even- 
indexed terms of D„, renumbered with the new index m {D {m = 3) = Dq, for instance). 
Comparing with expression ( |5^ above for the Zaheer-Zubairy states, we can see that the 
main difference in the present case is that the frequency of revivals is doubled: the A;*'' 

zz 



revival occurs around tf^ 



'^'^V'^^^'^'^ . ^ compared to = ^''^V^'^^+i ^y^q previous case. 
This difference is due essentially to the doubling of the separation between adjacent Rabi 
frequencies in the present case pTI], and is illustrated by comparing Figs. 2 and 4. 



V. CONCLUSION 

Intuitive pictures of the interaction between a two-level atom and an electric field com- 
monly involve the expectation that the atomic level populations must change as both sys- 
tems exchange excitations over the course of time |2^. This is due to the absence of further 



13 



atomic levels, which precludes the existence of destructive interference between different 
atomic transitions. However, in a fully-quantised interaction model such as the Jaynes- 
Cummings model, it is indeed possible to have states in which the atomic populations are 
completely or nearly completely trapped. This can be ultimately traced to the fact that 
the eigenstates of this model are entangled. We have shown that, by giving up the tradi- 
tional point of view based on the individuality of each subsystem and assuming instead one 
based on these entangled dressed states, it is possible to obtain a quantitative understand- 
ing of population trapping in this model. This is achieved via the introduction of a set of 
joint atom-field state variables, the 'weighted dressedness' distribution Dn (eq. |TUp. We 
have shown that, for general initial conditions, this distribution governs the evolution of the 
atomic inversion, assuming a role commonly attributed to the photon number distribution 
Pn (to which it reduces in particular cases). Using D„, we have obtained an upper bound 
to the amplitude of population oscillations that can be expected from a given state at any 
instant of its evolution. We have also been able to obtain an approximate analytical descrip- 
tion of the behaviour of the atomic inversion during partially suppressed revivals. We have 
found that in general the shape of revival envelopes is a direct reflection of the form of (a 
continuous interpolation of ) Dn- In the particular case of a field initially in a coherent state, 
this explains the appearance of a doublet structure in the revivals in the limit of greatest 
population trapping. 
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APPENDIX A: APPROXIMATE EXPRESSIONS FOR REVIVALS IN THE 
ATOMIC POPULATION INVERSION FOR GENERALISED INITIAL 

CONDITIONS 

In this Appendix, we calculate approximate expressions for the revivals in the atomic 
population inversion, which are valid for a wide variety of pure initial conditions of the 
atom-field system. These calculations generalise the method presented by Fleischhauer and 
Schleich [|], who assumed an initially factorised state of the form l^f)^ \i{^)p- The result 
they obtained for the k^^ revival is: 

where P„ is the photon distribution of state ■ In the case of more general initial con- 
ditions, including ones with atomic and atom-field coherence, we shall see that the form 
of the revival envelope depends not on the state's photon-number distribution, but on the 
'weighted dressedness' distribution D„ given in equation (|19|). 
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As was shown in section [II B 1| , the inversion, written in terms of dressed-state coordi- 
nates, has the form: 



(d-^ (t)) = -W_i + J2Dn COS (0„ (t)) 
n=0 



(A2) 



where 



0n (t) = <Pn (0) - (A3) 

This expression can be rewritten according to the Poisson Summation Formula (see 



Courant and Hilbert p2[, p.. 76 ) 



where 



Uk (t) 



dnD (n) cos (0„ (0) — fi„t) e 



2inkn 



To (t) = -D (0) COS (00 (0) - 2Xt) - w^i 



(A4) 

(A5a) 
(A5b) 



and where D (n) is any 'reasonable' function of a continuous variable n (continuous, differ- 
entiable, etc.), that interpolates between the values of Dn at the points where n integer. 

Noting that the sum in k in ( [^ ) extends to ±oo, so that the expression is invariant 
under k ^ —k, it is possible to substitute eq. (|A5a|) by: 



UJk (t) 



dnD (n) cos {(pn (0) -2Skin,t)) = (A6) 
dnD (n) cos (0„ (0)) cos {2Sk {n, t)) + 

oo 

dn sin (0„ (0)) sin {2Sk {n, t)) (A7) 

/•oo roo 

Re / dnDi (n) exp {2iSk {n, t)) + Im / dnD2 (n) exp {2iSk (n, t)) (A8) 

-'0 JO 



+ 



where we have defined 



Sk {n, t) = Tvkn — Xt^/n + 1 
Di(n) =D(n) cos (0,(0)) 
D2{n) =D(n)sin(0„(O)). 

In this way, the inversion may be rewritten (without any approximation) as 

oo 
k=—oo 



(A9a) 
(A9b) 
(A9c) 



(AlO) 



where 
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ul (t) = Re dnDi (n) exp {2iSk {n, t)) 
Jo 

POO 

ul (t) = Im / dnD2 (n) exp (2^5'^ (n, t)) 



(Alia) 
(Allb) 



Now, assuming the envelopes Di (n) and D2 {n) are sufficiently smooth if compared to the 
oscillating functions cos (25*^ (n, t)) , sin (25"^ (n, t)) , we may apply the method of stationary 
phases, approximating these expressions by: 



ool (t) ~ Di = Hk) Re 
ul (t) ^ 1)2 (n = Hk) Im 



where Uk is the point at which = : 





exp {2iSk {n - 


= nk)) X 




dnexp 


i 


{n - Uk f 




exp (2^5"^ (n 


= nk)) X 




dnexp 


i ^'^^ 


{n - nfc)^ 



(Al2a) 
(Al2b) 



2j.2 



nfc + 1 



(A13) 



We note that these expressions are invalid for A; = 0, since in this case the phase is always 
stationary (=0). Substituting in ( |749a| ) ,we obtain: 



Sk [n = Uk) 



d^Sk 



TV 



k 



2 = F 



(AM) 
(A15) 



Now, the integral in ( |y412a| , 6) may be written in terms of the Fresnel integrals p3 



C(x 



1 

IT Jo 



dy cos (x/"^^ dy 
^(^) ^ dy sin (y'^^dy 




C(x^oo) + C( \fFnk) +i{S 



For instance, for F > : 
dn exp {iF (■ 



[n - Uk) 

The asymptotic form of the Fresnel integrals for x — ^ 00 [^] is: 

C{x) 
Six) 



[X 



1 
2 
1 
2 



1 sin (x^) ^ / 1 

^ — - + O 

2ti X 

1 cos ix^ 



x^ 



27r X 



oil 



x" 



(A16a) 
(Al6b) 

00) + S[^nk)) 
(A17) 

(Al8a) 
(Al8b) 
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Assuming J\F\nk ^ 1 and taking these approximations to zeroth order, we have: 



dn exp {iF {n — 
Similarly, for F < 

j dn exp (j,F [n — 



[l + i) (vahd for F > 0) 



nk) 



TT 



2IFI 



(1-z) 



(A19) 



(A20) 



Using QAISI) and ( |A13| ) , condition J\F\nk ^ 1 becomes: 



At > 2 \k\ + J27r \k\ + 47r2A;2 



ou < 2 a/tt |A;| - J27r |A;| + Aix^k'^ 



(A21a) 
(A21b) 



for A; = 1, for example, this requires At ^ 17 

Substituting the asymptotic expressions (|A19|) , (|A20|) in ( |A12q| , h) , we obtain for /c > 
(or A; < 0): 



(t) ~Di(n = Hk) [cos (25fc|„^„J ^ sin (2,5^ 
it) -D2{n = Uk) [cos (25fcL=„J ± sin (^2Sk 

Thus, using (|A15|) and (|A10|) ,the inversion may be written: 



(A22a) 
(A22b) 



i ''0 


Di {n = Hk) 


V27rA;i / 


D2 {n = Uk) 



cos (2>Sfc|„^„J T sin {23^1^^^ 
) [cos (25fc|„^„J ± sin (2,Sfc|„^„^ 



+ 



(A23) 



(where the upper(lower) sign is valid for the terms with k > 0(< 0). Finally, substituting 
the value (^) for Sk\^ : 



it)) = E 



Xt 
2'nk^ 



k=—oo,k^O 
k=—ca,k^O 



^l + ^2L=„JCOS 



27rfc 



3 

27rA;2 



(27ck + 
sin (2vr/c + |^ 



± 



+ + + To 



\n=nk 



COS (0„ (0)) + sin (0„ (0))|_„J cos (fg ) ± 



± cos ((/)„(0))-sin((/)„ (0))| 



sm 



V 2-Kk ) 



(A24) 



Finally, since: 

(cos (x) + sin (x)) cos {y) ± (cos (x) — sin {x)) cos (y) = -\/2 cos ± y — — ^ (A25) 
then the approximate expression for the inversion is: 
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it)) 




D (n) cos (pn (0) ± 



27rA; 



TT 



?i+l= 



+ + CJg + To 



(A26) 



where the upper (lower) sign is vahd for the terms with A; > (< 0),and where tq and 
uJq + uJq are given by 



To (t) = -D (0) cos (00 (0) - 2Xt) - w.i 
ul + uJq (t) = uJo (t) = j dnD {n) cos ((pn (0) - 2\t\/n + l) 



(A27a) 
(A27b) 

The first of these two last terms represents the contribution to the inversion of the 
states in which the field is in a vacuum state. For the initial conditions which satisfy the 
approximations that have been made above, this term will normally be negligible (see below). 
The second term assumes non-negligible values only in the vicinity of t = 0. This is due 
to the fact that it is an integral over oscillating functions with different frequencies, which 
rapidly cancel each other out; also, since these frequencies form a continuum, they cannot 
re-phase substantially at subsequent times. Fleischhauer and Schleich have thus conjectured 
that zero-order terms of the Poisson Formula such as this always describe the first collapse 
of the inversion shortly after t = 0. 

The remaining terms ( A; 7^ ) each describe a modulated oscillation in the inversion, 
with an envelope given essentially by the shape of D (n) and assuming non-negligible values 

only in an interval around t ~ represents a value of n around the peak 

of D (n)). If D (n) is sufficiently narrow, the first few of these modulated oscillations will be 
well-separated in time, thus constituting an independent revival during which the inversion 
is described by 



it)) 



Xt 



v7rV2F 

When the initial state is of the form | (7)^ 



D (n) cos (j)n (0) ± 



TT 



IC 



2'Kk 4^ 
we have 



n+l= 



0„ (0) ^ TT , D (n) ^ P, 



71+1 



(A28) 



(A29) 



( Pn being the photon distribution of 
is recovered (eq. (2.86) in [|| ). 

This expression is valid as long as: 



so that Fleischhauer and Schleich's result ( |A1| ) 



a) The distributions D.„ cos (0„ (0)) and D„sin(0„(O)) vary slowly with n if compared 
with cosSk [n^t) = cos{nkn — \t\/n + 1). 



b) The value of t obeys conditions (|yl2lQ| , h) . This implies that, for the k^^ term of the 
sum above to describe well the k*^ revival, D (n) must assume its largest values in the 
region of n where: 



n 



1 > 



4 (Snk + iTT^k'^ + 2nk^/2 + ink 
47r2fc2 



4 + 



2Tvk \2 



V2 + Ank 



(A30) 
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Thus, the approximation can be expected to be good for all revivals if the initial state 
has at least 10 or so photons on average in the field. This will usually also ensure that 
the component tq, which depends on D (0), (eq. ( |A27aj )), can be ignored. 

It is straightforward to show that in the examples of section IV, where the dressed-state 
coordinates are given by expressions ( |5DD and where \a\ = 7, both of these conditions are 
satisfied. 



1. Variation for 'even-odd' states [where W2n~i = 0] 



In the case of 'even-odd'-type states such as {"^eo 7, 0)^ where only the dressed-state 
coordinates Wn with even index are non-null {w2n-i = 0), condition (a) above is violated 
and expression ( |A26| ) is thus invalid. Nevertheless, a similar analytical expression for the 
inversion can still be derived if eq. ( |^ ) is rewritten considering only the terms with even 
n. In this case, it is straightforward to show that, as long as the 'continuous versions' of 
distributions cos (02n (0)) and sin (02n (0)) are sufficiently smooth as a function of 
n, then the same stationary-phase method as was used above can be applied, resulting in: 




D (m) cos 



(A31) 



Here D (m) is a continuous interpolation of the even- indexed terms of Dn, renumbered 
with the new index m {D {m = 3) = Dq, for instance). Tq and + uJq are still given by 
( [/127a| , b) ,except one must substitute n ^ m. 

Once again, it is possible to show that the approximations realised in the course of 
obtaining this formula remain valid as long as D (m) assumes significant values only for 
2m ^ 10. (Thus, whenever the formula is applicable the term tq can be ignored). 
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FIGURE CAPTIONS 

FIG. 1. 'Dressed-state coordinates' for the Jaynes-Cummings model. Any state of 
the atom-field system can be represented in terms of the set of parameters 9n,4>n,Wn, Cn as 
defined in eqs. ([T0|) and (pT]). For each n, the first two of these can be pictured as forming a 
Bloch-type unit sphere, while the second two parametrise a circle of radius Wn- Under the 
effect of the Jaynes-Cummings Hamiltonian, w„ and 9n are constants of the motion, while 
(pn and undergo periodic motion at a frequency determined by the Rabi frequency fin- 

FIG. 2. Evolution of the atomic population inversion for an initial state where the field is 
in a coherent state and the atom is in an equally- weighed coherent superposition of |e) and 
\g), with relative phase ^ (see eq. (^)). The amplitude of the coherent state is a = |a| e^'^"', 
where |a| = 7 (49 photons on average in the field). In (a), the relative phase — 
equal to |, and the evolution follows the familiar collapse-revival pattern. As this phase 
difference is lowered to ^ (b) and finally zero (c), the revivals are quenched, and the atomic 
populations become effectively trapped (notice the scale change in (c) ). Simultaneously to 
this flattening, the revivals also assume a doublet structure. 

FIG. 3. 'Weighted dressedness' distributions D„ = w'^sinOn for states where the field is 
in a coherent state |a) with amplitude a = Te*"^" and the atom is in a coherent superposition 
COS7 |e) + e~*^ sin7 \g). In (a), 7 is fixed at | (equal weights), while the relative phase {ua—O 
varies from | (I) to j (II) to zero (III). The resulting flattening of D„ is due to the increasing 
dressedness of the components with the greatest weights Wn- In (b), (t'a — is fixed at zero, 
while 7 varies from | (1) to | (2) to ^ (HI)- In this case, the flattening corresponds to the 
matching of the points TT-max of maximum weight and n^i^ of maximum dressedness (see eqs. 



30| ). The limit of maximum flattening corresponds to the greatest degree of population 
trapping in the time evolution. The appearance of a doublet structure in this limit is due 
to the most important component of the state becoming completely dressed (sin6'„^^^ — >• 0). 
This is reflected in the shape of the quenched revivals, as can be seen in Fig. 2. 

FIG. 4. Evolution of the atomic population inversion when the atom-field system is 
initially in the entangled 'even-odd' state l^'^jo (a, 7, ,^)) given in expression (p^). The 
parameters in (a)-(c) are the same as those in Fig. 2. As in that quenching of 

revivals and appearance of a doublet structure is observed as (t'a — ~^ 0- The main 
difference in the present case is the halving of the interval between adjacent revivals, due to 
the doubled spacing between the Rabi frequencies present. 

FIG. 5. First and second revivals in the inversion for an initial state where the atom 
is in a coherent superposition of |e) and \g) and the field is in a coherent state \a) with 
amplitude a = 7. The top two graphs are close-ups of Figs. 2(a) and 2(c) and depict 
the numerical calculation of the Jaynes-Cummings sum (|18|) . The bottom two plot the 
approximate analytical expression (pOf), obtained from the Poisson Summation Formula. 
Despite a little distortion, agreement is seen to be good. Note that, at the tail end of (al) 
and (a2), the beginning of the third revival can already be seen interfering with the second 
one, while in (bl) and (b2) only the analytical expressions for first two revivals are plotted. 
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Fig. 2 




Fig. 3 




— 1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I 

20 40 60 80 100 



20 



(b) 




0.045 






0.040 






0.035 












0.030 






0.025 






0.020 






0.015 


1 * • 

/ • \ 




0.010 






0.005 
0.000 
-0.005 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 



100 



23 




It 




24 



